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Abstract 

Using in-situ transmission electron microscopy, we have directly ob¬ 
served nano-scale defects formed in ultra-high purity tungsten by low-dose 
high energy self-ion irradiation at 30K. At cryogenic temperature lattice 
defects have reduced mobility, so these microscope observations offer a 
window on the initial, primary damage caused by individual collision cas¬ 
cade events. Electron microscope images provide direct evidence for a 
power-law size distribution of nano-scale defects formed in high-energy 
cascades, with an upper size limit independent of the incident ion energy, 
as predicted by Sand et al. [Eur. Phys. Lett., 103:46003, (2013)]. Fur¬ 
thermore, the analysis of pair distribution functions of defects observed in 
the micrographs shows significant intra-cascade spatial correlations con¬ 
sistent with strong elastic interaction between the defects. 
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1 Introduction 


Low dose neutron irradiation and plasma exposure gives rise to major changes in 
mechanical and thermal properties of plasma-facing components for ITER and 
future generation fusion devices |T] . Tungsten, chosen as a plasma-facing mate¬ 
rial for the ITER divertor thanks to its high-temperature strength and thermal 
conductivity is expected to be significantly affected by neutrons produced 
from D-T fusion. In addition to transmutations @ng, neutrons initiate colli¬ 
sion cascades in which nano-scale defects are produced 0- Collision cascades 
provide a source that drives microstructural evolution of materials under irra¬ 
diation, resulting in the formation of complex dislocation and defect networks. 
The response of irradiated tungsten to mechanical deformation or temperature 
gradients differs significantly from that of the unirradiated material. 

Owing to the difficulty of producing intense fluxes of neutrons with a fu¬ 
sion energy spectrum, self-ion irradiation is used to simulate neutron damage, 
and there has been recent work investigating individual ion cascade events in 
tungsten both experimentally^ E] and theoretically^ [TO] [TTj J5] . As any ion 
moving with over the threshold displacement energy can stochastically set in 
motion other ions, the structure of collision cascades is fractal in nature P'2]. 
Furthermore, at high energies cascades branch into subcascades [13] [14]. This 
fractal nature suggests that defect statistics may be governed by a universal 
power law similar to other systems exhibiting self-organized critical behavior 
m • Evidence for such power-law statistics has been found recently in molecular 
dynamics (MD) simulations of cascades in tungstenpTj. Experimental observa¬ 
tions at room temperature are consistent with such primary damage predicted 
by MD and evolved with Object Kinetic Monte Carlo (OKMC) simulations^. 

From the primary damage state at 10ns after a cascade event to an ob¬ 
servation seconds later, the timescale gap bridged by OKMC is wide. Even 
if physically justified and carefully calibrated, OKMC based on semi-empirical 
laws for defect interactions and mobilities cannot by itself validate the initial 
defect statistics. This problem is compounded in self-ion irradiation, where the 
evolution is affected by the escape of radiation induced defects to the surface: 
what is observed in in-situ TEM at room temperature is merely the remnant of 
immobile defects. 

There is currently no experimental way to make observations of individual 
damage events over the nano- to milli- second timescale, but it is possible to ob¬ 
serve a dramatically slowed microstructural evolution using self-ion irradiation 
combined with in-situ transmission electron microscopy (TEM) at cryogenic 
temperature. In this work we observe defects formed in 150keV and 400keV 
W+ ion cascades in ultra-high-purity tungsten foil at 30K. At this temperature 
isolated clusters of point defects, be they vacancy or interstitial in nature, are 
generally considered to be immobile [IS] . 

We show that the distribution of sizes of visible defects measured experimen¬ 
tally and predicted by MD simulations of cascades in foils does indeed follow a 
power law, with 150keV tungsten ions producing defect distributions with an 
exponent in good agreement with the foil MD simulations. We show further 
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that 400keV tungsten ions show a slightly lower exponent in agreement with 
bulk MD simulations |11|. indicating a deeper penetration and reduced surface 
effects. We also examine the defect pair correlation function and demonstrate 
that the visible defects are not distributed homogeneously in space, but rather 
show significant correlation between their positions. An estimate for the elastic 
energy of interaction between nano-scale dislocation loops of such correlations 
shows that for loops with diameter d > 4nm this energy typically exceeds 1 eV. 
Our observations also show there exists a maximum defect size, of order 700 
point defects, which we attribute to sub-cascade branching. 

2 Experimental observations 

Ultra-high-purity tungsten sheets (typically W > 99.996 wt%) supplied by 
Plansee Gruppe were used in this study. 3 mm discs punched from the sheets 
were mechanically thinned to around 100 /im, annealed in vacuum at 1673 K for 
20 hours and finally twin-jet electropolished for TEM in 0.5 wt% NaOH aqueous 
solution. The heat treatment produced a {100}<011> texture, and the aver¬ 
age grain size exceeded 10 /im. Grain boundaries and residual dislocations are 
therefore not expected to significantly influence radiation defect production or 
evolution. Irradiations were performed in situ on the IVEM-Tandem Facility 
at Argonne National Laboratory, with 150 keV and 400 keV W + ions at 30 
K respectively, up to ~ 10 16 W + m -2 at around 15° off the zone axis of (001) 
grains. Foils chosen for microscopy analysis were ~100 nm in thickness. 

Damage profiles calculated using SRIM2013 m monolayer damage calcu¬ 
lation with E d = 55.3eVjS] suggest that the peak damage for 150keV ions is 
0.014dpa at lOnm depth; for 400keV ions it is 0.018dpa at 16nm depth. This is 
below the fluence for cascade overlap [7], so features we see in the micrographs 
can be attributed to individual cascade events. Both ion energies are expected 
to produce cascades contained within the foil. Weak-beam dark-field conditions 
(g = 200, ~ 3.8g) were selected for imaging defects. Image simulations [l8j 
show that size of an image reflects the true defect size under imaging conditions 
characterized by a high excitation error, close to s = 0.2 nm -1 . At this defect 
size and with these imaging conditions, loops appear as single spots and do 
not show a ‘coffee-bean’ double-lobe image. Figures [T] and [2] show micrographs 
captured in situ and the final damage state. 

3 Simulation methods 

Molecular dynamics simulations of collision cascades resulting from self-ion irra¬ 
diation in W foils were performed with the PARCAS code [m Eol eh. Tungsten 
atoms with 150 keV kinetic energy were initiated above a (014) surface, 15° off 
the surface normal to replicate the experimental geometry and to avoid strong 
channeling directions, with the foil atoms initially at OK. A simulation box of 
48x48x65 nm was used, with periodic borders in the transverse directions and 
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Figure 1: A sequence of snapshots taken over Is irradiation with 400keV 
ions showing the space- and time- correlated appearance of nano-scale defects. 
Contrary to the case in iron, where the defects produced at low fluence are very 
small [0, individual cascade events in tungsten produce correlated clusters of 
defects up to 9nm diameter, which appear at full brightness in a single frame 
and then do not subsequently evolve. 
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Figure 2: Defect clusters produced in collision cascades of tungsten. Top: 
Weak-beam dark-field imaging (g = 200, ~ 3.8g) of defect clusters in (left, 
centre) 150 keV, 400 keV W + ion irradiated samples, up to a fluence of 
~10 16 W + m~ 2 at 30 K. Top right: an enlarged view showing spatially corre¬ 
lated defects. Bottom: computed generated interpretations of the position and 
size of the defects in the images above. Defects are coloured according to their 
maximum diameter, with the scale (in nm) given in the coloured bar. 
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Figure 3: Histogram of defect size (N) for both experiments and MD simu¬ 
lations. Loops smaller than 1.5nm were difficult to distinguish from the back¬ 
ground image intensity fluctuations; the smallest size experimental bin is the 
accumulated frequency for diameter d < 1.5nm. The dashed line shows the clus¬ 
ter power law found for 150 keV cascades in bulk tungsten HU. The histograms 
are normalized so that the area under the curves matches the total number of 
defects counted. 

open boundaries at top and bottom surfaces. Atoms along the periodic borders 
were thermostated to 0 K using a Berendsen thermostat [22]. A modified EAM 
potential [23] was used, with the repulsive part given by the universal ZBL po¬ 
tential HM- Electronic stopping was included for atoms with kinetic energies 
above 10 eV PH EU- Counting and identification of defects was performed by 
an automated Wigner-Seitz method HU. The simulation time (50ps) was suf¬ 
ficient for the final defects to be internally stable, and would require thermal 
migration for further recombination to occur. 


4 Results 

The size of the defects and their relative positions were extracted from exper¬ 
imental micrographs using an automated procedure detailed below. For easy 
comparison to MD simulation results, we assume that each spot is formed by 
a circular y 2 (111) Burgers vector dislocation loop with diameter d, containing 
a number of point defects N = ny/3 (d/(2ao)) , where clq = 0.3165nm is the 
lattice parameter of tungsten. Measured and simulated defect size distributions 
were collected as histograms with varying width bins (fig. [3|. 

Over the loop size range N = 50 — 500 point defects, corresponding to loop 
diameters 2-6nm in the micrographs, the observed defect size frequency from 
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Figure 4: Typical final state configurations of (left) 30keV and (right) 150keV 
collision cascades in tungsten foil simulated using MD. Red spheres denote va¬ 
cant lattice sites, blue spheres are self-interstitials and green are surface atoms. 


150 keV ions is fitted to a power law / ~ A/N s with A = 2.7 ± 1.4 and 
S = 1.8 ± 0.1. This is in good agreement with the distribution derived from 
150keV foil simulations, where A = 18 ± 7 and S = 1.85 ± 0.09. 400keV ions 
produce a defect size distribution with A = 3.4±2.3 and S = 1.6±0.2. This can 
be compared to the frequency distribution with A = 7.5±1.5 and S = 1.63±0.07 
found in MD simulations of cascades occurring in the bulk El- 

Lower energy, 30 keV foil irradiation simulations show that most impacts 
result in heat spikes at very shallow depths, causing liquid-like cascade cores 
to erupt through the surface. This produces considerable sputtering, as well as 
viscous flow of the liquid atoms onto the surface [58], resulting in the forma¬ 
tion of near-surface vacancy-rich zones. Such depleted zones are seen in field 
ion microscope observations [29], and are well reproduced by MD simulations 
[3D]. Such cascades produce only a few small interstitial defects, biassing the 
size distribution to smaller defects. The agreement between 400 keV experi¬ 
mental observations and bulk MD simulations suggests that surface effects are 
not predominant in 400 keV ion irradiation. Ions penetrating deeper give rise to 
heat spikes contained in the bulk of the material, resulting in damage which is 
statistically similar to bulk damage, with roughly equal numbers of interstitials 
and vacancies, and larger dislocation loops. Figure [4] shows the results of typical 
cascades from a 30 keV and 150 keV W ion. While only 5% of cascades from 30 
keV ions resulted in a heat spike contained beneath the surface, this proportion 
was 50% for 150keV ions. 

For loops containing more than N = 700 defects, there is a deviation from 
linearity on the log-log plot. This deviation is significant as it indicates an upper 
limit to the size of defects generated in the cascades. Such a limit must exist, 
as it is impossible for a finite energy incident ion to generate arbitrarily large 
dislocation loops. 
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No loops larger than lOnm in diameter (N = 1300) were observed in the 
150keV irradiation experiments, and only one loop over lOnm in diameter 
formed during 400keV irradiation. It is not unlikely that this single very large 
loop was generated by the coalescence of other loops. The upper limit of the 
size distribution is similar for both incident ion energies, and the deviation from 
a power-law occurs at the same point. Such a deviation in the large event size 
limit is not uncommon in self-organized critical phenomena, see for example Fig. 
8 of Ref. [15| . In our case, the likely reason for the low number of very large 
defects is the branching of cascades into sub-cascades. This branching occurs 
above a certain threshold cascade energy, which for tungsten is close to 150 keV 

mm- 

The 150keV ion irradiations produced on average 0.27 visible loops per in¬ 
cident ion. This can be directly compared to 0.03 visible loops per ion in UHP 
W at 300K[8], indicating that at room temperature 90% of loops are sufficiently 
mobile to be lost to the surface or annihilate each other. The average loop diam¬ 
eter at 30K is 2.6nm (N = 90), compared to 2.0nm at room temperature, owing 
to the preferential annihilation of large loops 0. The 400keV self-ion irradiation 
experiments produced 0.55 visible loops per incident ion, with average diameter 
d = 2.9nm, (N = 110). 

The number of visible (diameter >1.5nm) defects per cascade predicted by 
MD is greater than that seen in experiment; the foil simulations expected to 
produce 1.1 defects per ion if the power law is extended to N = 700. On the 
one hand, this is not surprising given the variability associated with the choice 
of interatomic potentials [21 EH. On the other hand, the short time scale of 
MD simulations means that even mobile defects have little time to migrate and 
annihilate, or reach the surface. The total number of defects found in a cas¬ 
cade simulation represents an upper bound to the possible number of surviving 
defects observed on experimental time scales. A further apparent feature of 
Fig-! is that the experimental distribution extends to significantly larger sizes 
than the distribution generated by MD. This unusual and interesting aspect of 
comparison between simulations and observations highlights the difference be¬ 
tween sampling in the two cases. We were able to analyse only 50 MD cascades, 
whereas the 150keV experiment comprised 10000 cascades, and the 400keV ex¬ 
periment 20000 cascades. In situ experiments therefore make it possible to 
observe events too infrequent to find in MD simulations. The average diam¬ 
eter of visible loops in the foil simulations is d = 2.6nm, in agreement with 
experiment. 

The radial distribution function gd(r), defined in eqn!as the probability of 
finding pairs of spots with diameters greater than d, separated by distance r is 
shown in fig. UK We see a maximum correlation at r = 4nm, for d = 2 — 4nm 
spots in the micrographs, which is the same as the scale (4nm) of the projected 
radial distribution of clusters with N >= 30 found in cascade simulations. Such 
large clusters occur only infrequently in MD simulations, making it difficult 
to draw a corresponding smooth curve for gd> o(r). The maximum of gd=o{f) 
for the MD simulations is 9nm, decreasing to 6nm if we exclude single point 
defects from the analysis. The average radius of gyration r gi also decreases 
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Figure 5: Analysis of spatial correlations of defect clusters obtained from mi¬ 
crographs of 400 keV W+ irradiated tungsten at 30 K. (a) The radial pair dis¬ 
tribution function gd{r) is plotted against the separation distance r for defects 
with diameters d > 0, 2,4, 6 nm. The solid lines are to guide the eye.(b) A char¬ 
acteristic elastic energy scale computed from the size and estimated separation 
of the loops (see text). 


from r g = lOnrn for N > 1 to r g = 4nm for N > 10. This shows that in 
simulations, larger defects form closer together than small defect clusters or 
point defects. 

From the relative positions of the loops in the image we can compute a typical 
elastic interaction energy. While it is not possible to judge the relative depth 
of loops, their orientations or Burgers vectors, we can nevertheless estimate the 
order of magnitude. The elastic energy of interaction between two defects at 
distances greater than the size of the defects is 

= (1) 

where R is the relative position of defects a and b, P-" 1 and Pj,j > and the 
defect dipole tensors, and Gjfc(R) is the elastic Green’s function. Simplified 
expressions derived from this equation are given in Refs. HUSH nu. For each 
loop in the 400keV irradiation images, the maximum value of this energy for each 
loop pair was computed. The result (fig. §>) illustrates that the elastic energy 
of interaction between pairs of loops is typically greater than leV for loops 
with diameters over 4nm. This agrees with the analysis of elasticity-mediated 
ordering of nano-scale dislocation loops EH EH, and is a direct experimental 
confirmation of the occurrence of elastic trapping of interacting dislocation loops 
in cascades found in OKMC simulations^?;. 
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5 Conclusions 


In this joint experimental-theoretical work we have studied characteristics of 
the primary stage of radiation damage in pure tungsten foils due to self-ions. 
By performing in situ irradiations at 30K, we have restricted the mobility of 
defects produced in the collision cascade. To a first approximation, we can claim 
all defects are generated in individual cascades and become immobilised rapidly 
after the cascade cools, and so what we see in the TEM is very close to the final 
stable configuration found in (an ideal) MD simulation. 

We have shown that defects generated in cascades can not be treated as 
isolated. Figure [5|i shows that the positions of dislocation loops are correlated, 
and fig. shows that this implies significant elastic forces between the loops. 

We have also shown that defect size distributions from our MD simulations 
are in very good agreement with experiment. In particular we have demon¬ 
strated that a power-law distribution of loop sizes holds for cascades in ion- 
irradiated tungsten foils as well as in bulk. The spatial distributions of simu¬ 
lated cascades show that large dislocation loops are generated in close proximity 
to each other, at distances corresponding to those seen in the TEM. 

These results show that the number of point defects produced per inci¬ 
dent ion is not the only relevant metric for the subsequent evolution of the 
microstructure. We have shown that the distribution of loop sizes and their rel¬ 
ative positions can be measured experimentally and are well-reproduced by MD 
simulation. This information is critical for determining whether elastic inter¬ 
actions are sufficient to bring defects together, forming an obstacle for further 
dislocation motion, or whether thermal fluctuation is sufficient for defects to 
find alternate sinks and be removed. 
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A Analysis of micrographs 

Defects visible in the micrographs were detected with a newly developed mod¬ 
ification of Ridler and Calvard’s isodata algorithm m , proven against an ex¬ 
tensive catalogue of measurements made by hand. Once obvious artifacts such 
as foil edges and screw dislocations are removed from the image, the algorithm 
proceeds in four stages: flattening, identifying, comparing, and counting. 

The flattening stage removes long-wavelength (10s to 100s nm) variations 
in image intensity due to variations in the foil thickness. For each pixel i in 
the micrograph with intensity fi, the average intensity ft in a region 50x50nm 
is computed, excluding bright pixels over the Ridler and Calvard threshold [31]: 
/thresh = 1 /2(/fore + /back)- Then the intensity of the pixel is adjusted to 

fi (fi - 1) ^ j. -/) ’ ( 2 ) 

where / is a target background intensity. This flattening is repeated, averaging 
over a smaller lOxlOnm region. 

Next all pixels with /* > f x = x /i (/thresh + /back) are identified. A spot 
is defined as a contiguous region within this subset which includes a single 
maximum with intensity fi > /thresh- 

A path may exist between two spots along which the intensity does not drop 
below f x . If this is identified, it is necessary to determine whether we indeed 
have two spots, or if it is merely intensity fluctuations. For each pixel on a 
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straight line connecting two maxima, if the intensity drops below 70% of the 
interpolated value, or below f x , the maxima are treated as two spots. If no 
point on the line drops below this threshold, the two maxima are merged. 

With the positions and intensities of all pixels in the spot, the maximum 
diameter can be computed as d = 4r+, where r+ is the square root of the larger 
eigenvalue of the tensor T = /j(fj - (f)) ® {n - (r))/ fi¬ 
ll this procedure is applied to a reference image of unirradiated sample, a 
number of spurious dim spots will be identified, corresponding to short- wave¬ 
length variations in intensity present in all TEM micrographs. The mean and 
standard deviation of spots in the reference image is stored as a function of spot 
size. A spot in an image is accepted if its intensity passes a % 2 test at the 5% 
significance level. 

With the position r a and diameter d a computed for each spot as above, 
counting is straightforward. We add to the analysis the pairwise distribution 
gd{r), defined as 

g d (r) = ^ H( - da ~ d ) H ( d b “ d ) s (\n - r a \ - r), (3) 

a,b^a 

i.e. a count of all pairs of spots which both have diameters greater than d and 
are separated by distance r. H(x) is the Heaviside function. The normalisation 
Z is chosen such that the expectation value of gd{r) is one if the spots are 
randomly placed, computed numerically by randomly repositioning the spots a 
large number of times. 
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